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Abstract. We present the first long-term optical and in- 
frared study of the optical counterpart to the source RX 
J0812. 4-3114, an X-ray pulsar with a Be type companion. 
During the period covered by the observations the profile 
of some Balmer lines changed from absorption to emission 
and back again to absorption. Contemporaneously, the in- 
frared magnitudes varied by more than 0.8 mag. This long- 
term variability is interpreted as the formation and subse- 
quent dissipation of the Be star's disc. The building up of 
the disc ended up in an active X-ray state characterised 
by regular outbursts occurring at 80 day intervals. The 
overall duration of the formation/dissipation of the disc is 
found to be <4.3 years. Optical spectroscopic and infrared 
photometric observations were used to refine the spectral 
type of the primary (B0.2IV) and to monitor the circum- 
stellar envelope around the Be star. UBVRI and uvbyfi 
photometric observations allowed the determination of the 
astrophysical parameters of the optical companion. 
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1. Introduction 

The star LS 992 is the optical counterpart to the X-ray 
source RX J0812. 4-3114. Such an association was origi- 
nally proposed by Motch et al. (1997), as a result of a 
systematic cross-correlation between the ROSAT All Sky 
Survey (Voges et al. 1999) and several OB star catalogues 
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in the SIMBAD database. The X-ray light curve of LS 
992/RX J0812.4-3114 is characterised by 31.88 second 
pulsations, while the X-ray spectrum is best represented 
by an absorbed power-law component with a exponen- 
tially cut-off (Reig & Roche 1999). In December 1997 the 
source made a transition from a quiescent state to a flaring 
state (Corbet & Peele 2000), in which regular flares sepa- 
rated by 80 day intervals were detected with the All-Sky 
Monitor (ASM) onboard the Rossi X-ray Timing Explorer. 
Corbet & Peele (2000) attributed the origin of these flares 
to the periastron passage of the neutron star, hence this 
periodicity was naturally associated with the orbital pe- 
riod. The peak X-ray luminosity during outbursts was es- 
timated to be of the order of 2.3 x 10 36 erg s _1 , assuming 
a distance to the source of 9 kpc. 

In the optical band, LS 992 appears as a relatively 
bright V ~ 12.4 B0-1III-V star (Reed 1990). The optical 
spectra revealed Ha in emission and evidence for short- 
term (~ weeks) V/R variability (Motch et al. 1997). The 
combination of an emission line B star plus a neutron star, 
whose existence is inferred from the characteristics of the 
X-ray emission, defines a Be/X-ray binary. Be/X-ray bi- 
naries comprise approximately 70% of the more general 
class of high mass X-ray binaries, the remaining ~ 30% 
contain evolved (luminosity class I and II) primaries. The 
high energy radiation is the result of accretion of the ma- 
terial expelled by the Be star onto the neutron star. 

In this work we present the first long-term study of the 
optical counterpart of RX J0812. 4-3114 and discuss the 
consequences of the correlated behaviour in the optical, 
infrared and X-ray wavelength bands. 

2. Observations 

Blue spectroscopy of the source was obtained on 31st Oc- 
tober 1999 using the 1.52-m telescope at La Silla Observa- 
tory, Chile. The telescope was equipped with the Boiler & 
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Table 1. Journal of the optical spectroscopic observations 



Date 


MJD 


Wavelength EW(Ha) 


EW(H/3) 






coverage (A) (A) 


(A) 


low-resolution spectra 


07-10-96 


50364 


3200-8000 -6.5±1.0 


-0.6±0.4 


05-02-98 


50850 


3800-7800 -16.8±0.8 


-1.4±0.5 


red-end spectra 


04-04-96 


50178 


6300-7000 +1.4±0.5 




03-02-98 


50848 


6300-7000 -17.0±1.0 




09-01-99 


51188 


6200-6900 -20.5±0.8 




blue-end spectra 


05-04-96 


50179 


4200-5000 


+1.6±0.4 


05-03-98 


50878 


3700-5200 


-1.4±0.3 


10-01-99 


51189 


3700-5500 


-2.2±0.3 


31-10-99 


51483 


3800-5000 


-1.8±0.1 


15-09-00 


51803 


3000-5300 


+2.1±0.1 



Chivens spectrograph + #32 holographic grating and the 
Loral 38 camera. This configuration gives a resolution of 
i~ 0.5 A/pixel. Measurements of arc line FWHM indicate 
a spectral resolution of ~ 1.4 A at ~ 4500 A. A second 
spectrum was obtained on 15th September 2000, this time 
using holographic grating #33, which gives a resolution of 
-1.0 A/pixel (FWHMk 3.0 A at ~ 4500 A). 

All other spectra were taken using the 1.9m telescope 
at SAAO, South Africa equipped with the Cassegrain 
spectrograph. The 1996 data were obtained using the Reti- 
con Photon Counting System (RPCS), whereas all subse- 
quent data were obtained using the SITe2 CCD detector. 
Details of the wavelength coverage achieved in each data 
set are given in Table 1. 

The optical photometric data were acquired at the 
SAAO 1.0- m telescope using the Tek8 CCD camera, em- 
ploying UBVRcIc: uvby(3 and Ha filters, at a plate scale 
of 0.35 arcsec pixel -1 for the January 1999 observations 
and 0.7 arcsec pixel -1 for the October 1996 observations. 
The CCD frames were pre-processed using IRAF's CCD- 
PROC package for bias subtraction, overscan removal and 
flat-fielding. Aperture photometry was performed in the 
STARLINK GAIA package, using observations of E-region 
standard stars (Menzies et al. 1989) for calibration to the 
Stromgren system. 

The infrared photometric data were obtained at two 
telescopes - the 1.9m equipped with the Mklll IR pho- 
tometer at SAAO, in South Africa and the 1.5m Tele- 
scopio Carlos Sanchez (TCS), Tenerife, Spain using the 
Continuously Variable Filter (CVF). 

3. Spectral type 

The spectra of LS 992 in the classification region are dis- 
played in Fig. [I]. The September 2000 spectrum shows 
no evidence for emission components, but it has a lower 
resolution and Signal-to-Noise Ratio (SNR) than the Oc- 



tober 1999 spectrum. It has been included, however, so 
that the obvious changes can be observed. In particular, 
we call attention to the very strong in-filling of the Hei 
lines at AA 4026, 4471, 4921 A typical of Be/X-ray bi- 
nary counterparts. In the October 1999 spectrum H/J is 
strongly in emission, with an asymmetric blue-dominated 
profile, while similarly asymmetric emission components 
are present in the cores of all Balmer lines down to He. 
Longwards of A ~ 4750 A the continuum is dominated by 
weak (presumably Fen) emission features. 

Photospheric lines are all shallow and broad indicating 
a fast rotator. However, given the high SNR of the spec- 
trum even the weakest metallic lines are detected. As is 
typical of Be/X-ray binary counterparts, the presence of 
weak Hen lines indicates a spectral class close to B0. In 
order to improve the classification we have compared the 
spectrum of the source with those of standard stars taken 
either at the same resolution with the same instrument or 
at higher resolution (from the list of Steele et al.1999). 

We have used the main sequence standards v Ori 
(B0V), r Sco (B0.2V) and o Per (B0.5V). The absence 
of He II A4200 A at this resolution indicates that the star 
is later than B0. If the star was in the main sequence, 
it would indicate a spectral type B0.5V or later. How- 
ever, the strength of Hen A4686 A when compared to 
the nearby On A4676 A feature and the detectability of 
Hell A4541 A argue against this later classification. This 
seems to indicate a higher luminosity class. 

In Fig. [l], the spectrum of LS 992 is shown together 
with those of the MK standard stars HD 48434 (B0III) and 
1 Cas (B0.5III). The smooth progression of all the Hen 
lines favours a spectral classification B0.2III for LS 992. 
However, some metallic lines seem weaker than expected 
and a luminosity class IV sounds more appropriate. Direct 
comparison with the spectrum of the B0.2IV standard ip 1 
Ori at a slightly lower resolution (from the digital atlas of 
Walborn & Fitzpatrick 1990) supports this interpretation. 
Therefore we will accept a spectral type B0.2IV, though 
B0.2III cannot be excluded. 



4. Astrophysical parameters 

The determination of the astrophysical parameters by 
means of photometric calibrations in a Be star is not as 
straightforward as in a non-emission line B-type star due 
to the presence of the surrounding envelope, which distorts 
the characteristic photometric spectrum. Thus one has to 
correct for both circumstellar and interstellar reddening 
before any calibration can be applied. We have made use 
of the iterative procedure developed by Fabregat & Tor- 
rejon (1998, hereafter FT98) to determine the intrinsic 
magnitudes. This method is based on the fact that there 
exists a correlation between the equivalent widths of the 
Balmer lines and the anomalies in the uvby photometric 
indices produced by continuum circumstellar emission. 
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optical observations 



Date 


MJD 


Telescope 


I 


Ierr 


R 


Rerr 


V 


Verr 


B 


Berr 


U 


Uerr 


05-10-1996 


50362 


S A AO 1.0m 


11.72 


0.02 


12.16 


0.02 


12.48 


0.02 


12.86 


0.02 


12.57 


0.09 


Date 


MJD 


Telescope 


y 


yerr 


b 


berr 


V 


verr 


u 


uerr 




,0err 


23-01-1999 


51202 


S A AO 1.0m 


12.52 


0.04 


12.89 


0.04 


13.11 


0.07 


13.33 


0.14 


2.47 


0.02 


infrared observations 


Date 


MJD 


Telescope 


J 


Jerr 


H 


Herr 


K 


Kerr 


L' 


L'err 






13-01-1996 


50096.56 


TCS 


11.39 


0.05 


11.14 


0.02 


10.84 


0.02 










06-02-1998 


50851.40 


S A AO 1.9m 


11.21 


0.01 


10.95 


0.01 


10.70 


0.02 


10.36 


0.21 






07-02-1998 


50852.51 


S A AO 1.9m 


11.25 


0.01 


11.01 


0.01 


10.66 


0.07 










08-02-1998 


50853.43 


S A AO 1.9m 


11.26 


0.02 


11.02 


0.01 


10.77 


0.02 


10.61 


0.63 






28-10-1998 


51115.76 


TCS 


11.67 


0.03 


11.31 


0.02 


11.03 


0.04 










05-01-1999 


51184.53 


S A AO 1.9m 


11.54 


0.02 


11.25 


0.07 


10.95 


0.02 










22-01-2000 


51566.59 


TCS 


12.01 


0.08 


11.89 


0.04 


11.43 


0.01 










22-01-2000 


51566.593 


TCS 


12.04 


0.33 


11.78 


0.14 


11.25 


0.02 












4300 
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Fig. 1. The spectrum of LS 992 in the classification region (middle) is compared to those of the MK standard stars 
HD 48434 (B0III, bottom) and 1 Cas (B0.5III, top). The inset shows the smooth progression of Hen AA 4541, 4686 
A (marked with arrows) when compared to the neighbouring Si in A4553 A and O n A4676 A. All spectra have been 
smoothed with a a — 0.8 Gaussian filter, normalised by division into a spline fit to the continuum and offset for display. 



4-1. Extinction, absolute magnitude and distance 

By applying the FT98 method we are able to disentan- 
gle the relative contributions of the interstellar and cir- 
cumstellar extinction: E is (b - y) = 0.433 ± 0.033 and 
E cs (b - y) = 0.051 ± 0.025, respectively. The errors re- 
flect the accuracy of the FT98 procedure. Using the rela- 
tion E 1S (B — V) = 1.35 E' s (b — y) for interstellar extinc- 
tion (Crawford & Mandwewala 1976) and E CS (B - V) = 
1.20 E cs (b — y) for circumstellar extinction (Fabregat & 
Reglero 1990) we find a total value for the extinction of 
E(B — V) = 0.65 ± 0.05, which is in good agreement with 



the values of 0.69 (from optical photometry) and 0.8±0.2 
(from interstellar spectral lines) reported by Motch et al. 
(1997). Correcting the B — V index for extinction we find 
-0.27, consistent with a B0.2IV star (Wcgner 1994). 

Although the (3 index (Crawford & Mander 1966) pro- 
vides a measure of the luminosity for O and B stars it 
is strongly affected by emission. Thus it cannot be used 
to calibrate the stellar luminosity. In view of the good 
quality of our spectrum we prefer to assume the lumi- 
nosity class derived from the spectral analysis and take 
My = —4.3 ± 0.5 as would correspond for a B0. 2 giant- 
subgiant star. The error has been arbitrarily taken so that 
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Table 3. Astrophysical parameters of LS 992 



Spectral type 
E(B-V) 

Toff 

Radius 

Mass 

M v 

Distance 

BC 

A/boi 



B0.2IV-III 
0.65±0.05 
28000±2000 K 
10±2 R e 
17±3 M Q 
-4.3±0.5 
8.8±4.0 kpc 
-2.8±0.2 
-7.1±0.5 
3.7±0.2 cm 2 s~ 
240±20 km s -1 



it includes disagreements between different My calibra- 
tions. The distance, estimated from the distance-modulus 
relation is found to be d ks 8.8 kpc. 

4-2. Effective temperature 

A B0.2IV star has a temperature close to 30000 K (Pana- 
gia 1973, Lamers 1981, Vacca et al. 1996). In O and B 
type stars the index Co is closely related to the effec- 
tive temperature (Crawford 1978). Using the calibration 
log T eB = 0.186cg - 0.580c + 4.402 (Reig et al. 1997a) 
we obtain T c g — 25400 ± 4500 K, where the error was 
obtained by propagating the photometric error. Due to 
the uncertainty in the u measurement, hence Co, the value 
of the temperature is not well constrained. Also, possi- 
ble residual circumstellar emission may result in slightly 
larger cq values. If the intrinsic colour derived above, 
(B — V)q = —0.27, is used in combination with the cali- 
bration of Gulati et al. (1989) a value of 26400 ± 2300 K 
is found. Whereas if Eq.(ll) of Napiwotzki et al. (1993) is 
applied, T e g k 29500 K. As a compromise, we adopted a 
temperature of 28000 ± 2000 K 

Directly related to the effective temperature there is 
the bolometric correction, BC = Mboi — My. The calibra- 
tions of BC in terms of T cS result in BC = -2.8 ± 0.2 
(Balona 1994, Vacca et al. 1996). 

4-3. Radius 



In order to estimate the radius of LS 992 we have re- 
written the luminosity equation L = Aira i? 2 T c 4 ff in terms 
of the visual brightness parameter Fy (Barnes & Evans 
1976) and the (b — y)o colour index. With the value of 
the reddening derived in Sect. 4.1 we obtain (b — y)o = 
—0.115 which gives Fy — 4.166 ± 0.026 according to 
Eq. (6) of Moon (1984), where the error is that of the 
least-squares fit of the calibration. Substituting in Fy — 
4.236 - O.lMy - 0.51og(i?/i? s ) (Moon 1985) we finally 
obtain R = 10.0 ±2.6 Rq. The final error was obtained by 
propagating the errors through the luminosity equation. 
This value compares well to the average radii in terms of 



MK spectral types calculated by Moon (1985), who give 
12 R Q and 7.6 R for B0III and B1III stars, respectively. 
And also to those of Vacca et al. (1996) who derived 8 and 
14.8 R for B0.5V and B0.5III stars, respectively. 

4- 4- Mass and Gravity 

With the values of the temperature and bolometric lu- 
minosity, we can estimate the mass from evolutionary 
models. We have used the interpolation formula given by 
Balona (1994) for the Claret & Gimenez (1992) models 
and found M = 17±3 M Q . Once we know the mass and 
radius of the star we can calculate its gravity by means of 
log g = 4.44 + log M - 21ogi?, where M and R must be 
given in solar units. The result is log<? = 3.7 ±0.2 cm s _1 . 

4-5. Rotational velocity 

The projected rotational velocity v sini can be estimated 
from the width of the spectral lines. However, in practice, 
it is difficult to find clean spectral lines, i.e., not affected by 
circumstellar emission. The spectrum taken in September 
2000 shows Balmer lines in absorption, indicating a disc- 
loss phase. In addition, He I lines are mainly formed in the 
star's photosphere and are expected to have the smallest 
disc contribution. We have used the He I A4026, He I A4144 
and He I A4388 lines and obtained FWHM= 6.1,5.9 and 
7.0 A, respectively. 

From Buscombe's (1969) approximation 



FWHM 



c 2A (ln2) 1 /2 



(1) 



where Ao is the rest wavelength of the line and FWHM 
represents the full width at half maximum corrected for 
instrumental broadening, we obtain v sinz =240 km s _1 , 
220 km s" 1 and 260 km s" 1 for He I A4026, A4144 and 
A4388 respectively. Similar values are found using the fits 
of Steele et al. (1999). 



5. Long-term variability 

Prior to 1998 LS 992/RX J0812.4-3114 was in an X-ray 
dormant state as indicated by the fact that the ASM flux 
was consistent with the X-ray background flux, ~ 0.1 c/s 
(Fig |). At the end of 1997 (MJD ~ 50800) the source un- 
derwent a transition from this inactive state to a flaring 
state, consisting of regular outbursts occurring at approx- 
imately 80 day intervals (Corbet & Peele 2000). At the 
onset of the X-ray activity the optical companion was in 
a bright infrared and optical state. 

In Be/X-ray binary systems the optical primary has a 
circumstellar disc, which is presumably made up of mate- 
rial expelled by the B star. A neutron star orbits around 
the Be star in a moderately eccentric orbit. The ultimate 
physical process which leads to mass ejection in Be stars 
is unknown but it is probably one (or a combination) of 
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Fig. 2. IR magnitudes, H/3 equivalent width and ASM Fig. 3. The large line profile variability indicates major 
light curve for the period April 1996-September 2000. The structural changes in the disc of the Be star. The contin- 
activity of the source increased suddenly at around 1997 uum has been normalised to unity and the spectra offset 
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these: radiation pressure, centrifugal force, thermal ex- 
pansion and hydromagnetic waves (Strafella et al. 1998; 
Balona 2000). This circumstellar disc is responsible for the 
two defining observational properties of a Be star, namely 
the infrared excess and the emission lines. The mechanism 
that powers the X-ray source is accretion onto a neutron 
star and the fuel is the material in the disc. 

In this scenario we would expect some type of correla- 
tion between the X-ray and the optical/infrared emission. 
Fig H shows the ASM and infrared light curves together 
with the evolution of the H/3 equivalent width. When the 
circumstellar disc becomes large enough the neutron star 
is able to interact with it and produce X-rays. As can be 
seen in Fig || the equivalent width of the H/3 line, EW(H/3) 
increased slowly during the years prior to the onset of 
the X-ray activity, indicating the building-up of the disc. 
Spectral hydrogen lines were in absorption in April 1996, 
in emission in October 1996 and back to absorption in 
Semptember 2000. The absorption profile indicates that 
the Be star's disc was at some point non existent. Subse- 
quently, the disc began to grow until it was large enough 
for the neutron star to interact. The source was X-ray ac- 
tive in December 1997. Thus, the build-up of the disc took 
< 18 months. The source reached an optical and infrared 
maximum and then began to decrease. 



The issue of how far/big the disc can grow was dis- 
cussed by Reig et al. (1997b), who concluded that the 
disc may be truncated by the continuous passage of the 
neutron star once the disc radius is comparable to the pe- 
riastron distance. This effect is expected to be stronger in 
systems with short orbital periods. Assuming that the 80 
day periodicity found in the X-rays (Corbet & Peele 2000) 
represents the orbital period, the diagram EW(Ha)-P or .f, 
(Reig et al. 1997b) predicts EW mol ~ 20-25 A, i.e, close 
to the maximum value reported in our long-term monitor- 
ing of the system (Table |]J). Thus, it seems reasonable to 
think that the disc in LS 992/RX J0812.4-3114 reached its 
maximum size by the time that the X-ray emission began. 

The onset of the X-ray activity marked the beginning 
of the disc-loss phase. This result can be understood since 
the material from the disc constitutes the fuel that pow- 
ers the X-ray source. The disc is also responsible for the 
infrared excess observed in Be stars, which is attributed 
to recombination radiation from ionized hydrogen (free- 
free and free-bound processes). Any major change in the 
properties of the disc will translate into the infrared mag- 
nitudes. Soon after the onset of the X-ray emission, the 
infrared magnitudes began to decrease and by the time H/3 
line went back into absorption J, H and K had decreased 
by more than 0.8 magnitudes. This relationship between 
X-ray and infrared emission also demonstrates that the 
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compact object exerts a notable influence on the fate of 
the Be star's disc. 

Overall, the entire episode of formation and dissipation 
of the disc, considering it as the time elapsed between two 
consecutive absorption line phases, took < 4.3 years (Fig 
||). It is interesting to note that the Ha line was showing 
little emission (EW(a)=-4.6 A) back in April 18, 1992 
(Motch et al. 1997), that is ~ 4 years before the first 
measurement reported here. That state must correspond 
with the first stages of a previous active state, and implies 
that the 4.3 yr period is repeatable. In other words, the 
formation and loss of the disc represents a cyclic event in 
LS 992/RX J0812.4-3114. 

The X-ray periodic flaring behaviour is common 
among Be/X-ray binaries and the flares, which normally 
peak at around periastron, are referred to as type I out- 
bursts. They differentiate from type II outbursts which 
are major increases (10 2 — 10 3 ) in X-ray flux and tend to 
peak at later orbital phase. Many of the giant outbursts 
are detected in the middle or at the beginning of a se- 
quence of normal outbursts (Bildsten et al. 1997). The 
case of LS 992/RX J0812.4-3114 is different since a se- 
ries of normal outbursts followed immediately the inactive 
phase. In this respect LS 992/RX J0812.4-3114 resembles 
EXO 2030 +375. This 46 d orbital period Be/X-ray binary 
was in an active X-ray state, between April 1991 (when 
BATSE became operative) until August 1993. It switched 
off for three years to reappear in April 1996, showing peri- 
odic X-ray outbursts at 46 d intervals (Reig & Coe 1998a; 
1998b). The similar X-ray luminosities of the flares (2- 
4 x 10 36 erg s _1 ) stresses the resemblance of these two 
systems. 

While the correlation between optical and infrared 
emission is the expected behaviour for a growing disc the 
long-term monitoring of Be/X-ray binaries is providing 
new insights into the way the neutron star affects the 
properties of the disc and the generation of high energy 
radiation. The traditional idea of expanding discs is be- 
ing replaced by disc dissipation/reformation cycles which 
control the X-ray behaviour. Once the disc reaches a cer- 
tain size, which presumably is limited by the continuous 
passages of the neutron star, the disc becomes more dense 
and unstable and material is fed onto the neutron star. 
Negueruela et al. (2000) have pointed out that disc trun- 
cation by the neutron star may indirectly cause the warp- 
ing/tilting of the Be star's disc. They applied the viscous 
decretion disc model (Negueruela & Okazaki 2000) to the 
Bc/X-ray transient 4U 0115+63/V635 Cas to explain its 
X-ray and optical/IR behaviour. 4U 0115+63/V635 Cas 
presented a sizable disc ~ 6 months after being observed 
discless and was X-ray active only 13 months after a pho- 
tometric low. These timescales compare very well with 
those in LS 992/RX J0812.4-3114. 

The asymmetry of the H/3 line profiles in LS 992 and 
the variability of its blue (V) and red (R) peaks (Fig ||) are 
reminiscent of the so-called V/R variability, and point to- 



ward a perturbed density distribution in the disc. The per- 
turbation is associated with the existence of global oscil- 
lation modes propagating though the disc (Okazaki 1991, 
1997; Papaloizou et al. 1992). The V/R pattern developed 
shortly after the disc build-up began, supporting the idea 
that the discs are intrinsically unstable against density 
perturbations. 

6. Conclusion 

We have witnessed the appearance and disappearance of 
the circumstellar envelope in the Be/X ray binary LS 
992/RX J0812.4-3114. Using UBVRI and uvby/3 pho- 
tometric observations we have derived the astrophysical 
parameters of the Be/X-ray binary LS 992/RX J0812.4- 
3114. The optical companion turns out to be a B0.2IVe 
star showing marked optical and infrared variability on 
timescales of months to years. The reversal from positive 
to negative and then back to positive of the equivalent 
width of the H/3 line, the large amplitude of variation of 
both, the Ha equivalent width (AEW(Ha) w 20 A) and 
the infrared magnitudes (AJ « AH « AK ss 0.8 mag) 
are interpreted as the formation and dissipation of the 
Be star's disc. The long-term monitoring of LS 992/RX 
J0812. 4-3114 constrains the timescale for formation and 
dissipation of the disc to about 4.3 years. Measurements 
found in the literature seem to indicate that this episode 
is cyclic and basically supports the picture of dynamically 
unstable discs, that is, the discs go through cycles of dissi- 
pation/reformation which control the X-ray behaviour of 
the source. 
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